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1 .  Introduction 

This  paper  concerns  Hamiltonian  and  non-Hamiltonian  perturbations  of 
integrable  two  degree  of  freedom  Hamiltonian  systems  which  contain  homo¬ 
clinic  and  periodic  orbits.  Our  main  example  concerns  perturbations  of  the 
uncoupled  system  consisting  of  the  simple  pendulum  and  the  harmonic  oscilla¬ 
tor.  We  show  that  small  coupling  perturbations  with,  possibly,  the  addition 
of  positive  and  negative  damping  breaks  the  integrability  by  introducing 
horseshoes  into  the  dynamics. 

We  begin  with  an  unperturbed  n  +  1  degree  of  freedom  Hamiltonian  in 
canonical  coordinates  q  =  (q\  ....  qn),  p  =  Cp-j ,  •••»  Pn)»  x,  y  of  the 
form 

H°(q,p,x,y)  =  F ( q , p )  +  G(x,y)  (1.1) 

Starting  in  §3,  we  will  assume  n  =  1 ,  but  for  some  of  the  development  n 
can  be  arbitrary.  Allowing  x  and  y  to  be  multidimensional  will  be  the 
subject  of  another  publication. 

We  shall  assume  that  G  admits  action-angle  variables;  i.e.  there  is 
a  canonical  change  of  coordinates  to  (6,1)  such  that  9  is  2it  periodic, 

I  >_  0  and  G  becomes  a  function  of  I  alone;  we  write  G(I)  for  this 
function  and  assume  that 

G(0)  =  0,  fl(I)  =  G ' ( I )  >0  for  I  >  0  (1.2) 

Note  that  (1.2)  implies  the  existence  of  G-1. 

The  equations  of  motion  are 


n1  -  3F 

q  '  3Pi  ’ 

•  3  F  •  % 

p,  =  -  i  =  1,  ...,  n 

1  Sq1 

(1.3) 

e  =  fid). 

i  =  o 

(1.4) 

3 

We  shall  assume  that  the  system  (1.3)  contains  a  homoclinic  orbit  (q(t  -  tQ), 
p(t  -  t  ))  joining  a  saddle  point  (qg,  Pg)  to  itself.  Of  course 

(1.4)  contains  the  2Tr-periodic  orbits  6(t)  =  0g  +  tft(Ig),  I(t)  *  IQ. 

Thus,  for  the  system  (1.3)-(1.4),  we  have  orbits  which  are  the  products 
of  the  homoclinic  orbits  and  the  periodic  orbits.  (See  Figure  1.)  [The 
case  in  which  F  has  a  heteroclinic  orbit  may  be  treated  by  similar  methods.] 


(1.7) 


<t>(t)  =  ±2  arctan[sinh(t  -  tg ) 
v(t)  =  ±2  sech(t  -  tg)« 

We  deal  with  Hamiltonian  perturbations  of  (1.1)  in  Sections  2,  3  and 
4.  We  assume  that  our  perturbed  Hamiltonian  depends  on  a  small  parameter 
e  in  the  form 

He(q,P,8,l)=  F(q,p)  +  G(I)  +  eH1 (q ,p,0,I) ,  (1.8) 


where  H1  is  smooth  and  2tt  periodic  in  e.  We  shall  show  that  a  Poincare 
map  associated  with  He  contains  Smale  horseshoes  on  each  energy  surface 
for  e  small  and  H1  satisfying  certain  conditions. 

The  equations  of  motion  corresponding  to  He  are 


JL  +  C^ 

9p.  e3p. 


c  idi. 
T  E  i  * 
3q  aq1 


0  »  J2(I)  +  e 


=  1, 


(1.9) 


Our  method  for  finding  horseshoes  involves  the  Melnikov  function  tech¬ 
nique  that  has  been  used  in  Melnikov  [1963],  Arnold  [1964],  Holmes  [1979, 
1980],  Holmes  and  Marsden  [1981]  and  Greenspan  and  Holmes  [1981]  to  show 
the  existence  of  transverse  intersections  of  stable  and  unstable  manifolds 
and  hence  the  existence  of  horseshoes.  The  Melnikov  technique  is  used  after 
the  system  has  been  reduced  to  a  non-autonomous  single  degree  of  freedom 
system  (as  in  Whittaker  [1959]  Ch.  12,  and  Birkhoff  [1966],  Ch.  VI,  §3). 

In  particular,  in  Section  4  we  prove  that  the  pendulum-oscillator  (1.6) 
develops  a  horseshoe  on  each  energy  surface  near  the  value  H  =  1,  when 
it  is  perturbed  using  the  coupling  term 


(1.10) 


H^(4>,v,x,y)  =  j(x  -4>)2 

Churchill  [1980]  suggested  the  possibility  of  this  approach  but  did  not 
examine  any  specific  examples.  Section  5  concerns  the  more  delicate  case 
in  which  (1.4)  is  given  an  additional  non-Hamiltonian  perturbation.  We 
prove  that  at  least  one  of  the  horseshoes  persists  under  this  perturba¬ 
tion  provided  there  is  a  suitable  energy  transfer  mechanism.  In  Section  6 
we  apply  this  theory  to  the  pendulum  oscillator  example  once  more. 

In  another  paper  (Holmes  and  Marsden  [1981b])  we  use  these  methods 
to  address  the  question  of  nearly  integrable  multi  degree  of  freedom  systems 
and  Arnold  diffusion  (cf.  Arnold  [1964]).  Holmes  and  Marsden  [1981c]  treats 
Hamiltonian  systems  with  symmetry  in  which  (part  of)  the  phase  space  is 
the  coad joint  orbit  of  a  Lie  group.  This  provides  a  natural  framework  in 
which  to  consider  non-integrable  pertrubations  of  rigid  bodies. 

In  many  examples  of  physical  interest,  such  as  weakly  nonlinear  prob¬ 
lems,  the  unperturbed  system  H  *  F(p,q)  +  G(I)  does  not  possess  a  homo¬ 
clinic  orbit,  but  some  averaged  system,  after  truncation,  does  have  homo¬ 
clinic  orbits  (cf  McGeehee  and  Meyer  [1974]).  In  such  cases  the  Melnikov 
function,  computed  with  the  use  of  second  order  terms  normally  neglected 
in  averaging,  is  typically  exponentially  small  and  conclusions  on  the  in¬ 
tersections  of  manifolds  do  not  immediately  follow  without  a  careful  study 
of  the  errors.  The  elastic  pendulum  in  the  limit  of  a  very  stiff  rod,  with 
linearized  frequency  and  Hamiltonian 

v2  ITT  fe) 

H  =  -y  -  cos  <j>  +  ml  -cj—  sin  -  cos  <p,  (1.11) 

also  falls  into  this  class.  The  study  of  such  systems  is  planned 
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for  a  future  publication.  The  problems  of  the 

motion  of  four  point  vortices  treated  by  Ziglin  [1980]  and  the  motion  of 
charged  particles  in  the  earth's  magnetic  field  (see  Braun  [1981])  possess 
related  difficulties. 

We  expect  that  the  methods  developed  here  will  be  applicable  to  a  num¬ 
ber  of  Hamiltonian  systems  exhibiting  complex  dynamics.  Two  examples  that 
seem  to  involve  homoclinic  phenomena  are  the  Henon-Heiles  system  (see 
Churchill,  Pecelli  and  Rod  [1979])  and  the  mixmaster  model  in  cosmology 
(Barrow  [1981]).  The  results  in  §5  should  also  enable  one  to  deal  with 
nearby  systems  with  forcing  and  dissipative  terms. 

For  other  papers  in  which  horseshoes  are  found  in  two  dimensional 
mappings  by  very  different  techniques,  see  Devaney  and  Nitecki  [1979]  and 
Tresser,  Coullet  and  Arnoedo  [1979]. 

Acknowledgements .  A  number  of  helpful  comments  were  kindly  supplied 
by  Allan  Kaufman,  David  Rod,  and  Alan  Weinstein. 
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2.  The  Reduction  Method 

We  now  recall  how  to  reduce  the  n  +1  degree  of  freedom  system  (1.9) 
to  an  n  degree  of  freedom  non-autonomous  system.  This  is  a  special  case 
of  the  general  reduction  procedure  by  which  a  Hamiltonian  system  with 
syinnetry  is  reduced  to  another  Hamiltonian  system  with  fewer  degrees  of 
freedom.  The  standard  reference  is  Whittaker  [1959,  Ch.  12];  see  also 
Birkhoff  [1966]  and  Churchill  [1980].  The  case  of  concern  in  this  paper 
is  the  symmetry  of  time  translations, with  energy  being  the  corresponding 
conserved  quantity.  The  procedure  is  also  a  special  case  of  that  of  Marsden 
and  Weinstein  [1974]  in  the  context  of  time-dependent  mechanics,  as  in 
Abraham  and  Marsden  [1978,  §5.1], 

Energy  is  conserved  for  (1.9),  so  we  consider  the  equation 

He(q,p,e,I)  =  h  .  (2.1) 

gHe 

Now  =  n(I)  +  e  -gj-  .  On  any  compact  subset  of  (q,p,0,i)  space  not 

3He 

containing  1=0,  we  can  choose  e  small  enough  so  that  -gj-  >  0,  since 
,:(I)  >  0  for  1^0,  by  assumption.  Thus,  in  such  a  region,  we  can  solve 
(2.1)  for  I  to  obtain 

I  =  Le(q,p,6,h)  *  (2.2) 

Now  define  and  iJ  by  writing 

Le(q,p,9,h)  =  L°(q,p,h)  +  eL1 (q,p,0,h) +  0(e2)  ,  (2.3) 

2 . 1  Proposition.  We  have 

L°(q,p,h)  =  G_1(h  -  F(q,p))  y  (2.4) 
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and 


1 


L  (q,p,0,h) 


.  HUq.P.Q,  L°(q ,pth) )  , 


TJ 


n(L  (q.p.h)) 


(2.5) 


Proof.  Substituting  (2.2)  into  (2.1)  gives 

F(q,p)  +  G(L°  +  eL1  +  0(e2))  +  eH^q.p.e.L0  +  eL1  +  0(e2))=  h 


i.e.  [F(q,p)  +  G(L°)  -  h]  +  J2(L°)(eL1)+  eH1  (q,p,0,L°)  =  0(e2) 


The  and  terms  of  this  expression  give  (2.4)  and  (2.5). 


Having  passed  to  the  level  set  H  =  h  and  thereby  eliminated  I,  we 
now  eliminate  thevariable  conjugate  to  H,  namely  t.  (In  reduction,  an 
even  number  of  variables  is  always  eliminated).  Since  12(1)  >  0  and  He 
is  not  explicitly  t-dependent,  and  6  is  (for  small  e)  an  increasing  func¬ 
tion  of  t,  we  can  eliminate  t  by  inverting  e  =  0(t)  and  expressing 

i 1  i 

q  and  p  as  functions  of  0.  We  write  q  =  dq  / de  and  pj  =  dp./de, 

i  =  1 ,  . . . ,  n  -  1  so  that 


f '  g1  3He  /  3He 

q  *  t  -  3fr/  -31  ' 


and 


Pi 


.  p'  _  3H£  Id He  , 


(2.6) 


aq 


However,  implicit  differentiation  of  (2.1)  gives 


3He  .  3He  3Le  _  n 
Sq1  31  Sq1  * 


3H£  +  3H^  31/  q 
31  3pi  ‘  ' 


(2.7) 


9pi 


I 
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I 


Substituting  (2.7)  into  (2.6)  we  get 

i '  3Le 


3Pi  ' 


and 


j!  =  C 

1  W 


Using  (2.3),  this  becomes 


i ' 

q 


3L°  SL1  .  2, 

¥7  '  E  3PT  °U  )s 


4*  t^*0(e2)- 

3q'  Sq1 


Since  L°  depends  only  on  (q,p),  but  L1  depends  on 
system  (2.9)  has  the  form  of  a  2iT-periodical ly  perturbed 
freedom  Hamiltonian  system.  For  n  =  1,  (2.9)  becomes  a 
equation.  This  is  exactly  the  situation  which  occurs  in 
tor  problem. 


L  * 


(2.8) 


(2.9) 

q,p  and  0,  the 
n  degree  of 
forced  oscillator 
our  pendulum-oscilla 
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3.  Melnikov's  Method:  The  Existence  of  Horseshoes 

For  n  =  1 ,  the  system  (2.9)  is  in  the  form  analyzed  by  Melnikov 
[1963],  Holmes  [1979,  1980]  and  Greenspan  and  Holmes  [1981].  [For 
vectorial  or  n  >  2,  analogous  technique  were  developed  by  Arnold  [1964] 
and  Holmes  and  Marsden  [1981b]  and  will  be  the  concern  of  a  subsequent 
paper.] 

For  e  =  0,  the  system  (2.9)  reduces  to  (1.2)  and  thus  also  contains 
a  homoclinic  orbit.  The  Melnikov  method  involves  integration  of  the 
Poisson  bracket  around  the  homoclinic  orbit  of  the  unperturbed 

system.  Let  us  first  use  (2.4)  and  (2.5)  to  express  (L^L1  }  in  terms 
of  {F.H1}. 

3.1  Proposition.  Holding  0  and  h  fixed,  we  have 

aV)  -  — Vt  <m'>  (3.i ) 

fea°n2 


Proof.  Using  (2.4),  we  have 


3q 


9F 

“  l 

■ i 


l  3q 


anA  3LU  _  3F] 

and  apT  -  (G  >  *  55T 


while  (2.5)  gives 


3L '  = 
3q' 


na^y 


3H1  t  3H1  3L° 

(Sq1  31  3q1-J 


1  i.l 3L° 

- — x  H  fl  — r- 

[na°)r  aq1 


and 


3L 
3  Pi 


^H1  3H1  3L^ 

3p1  "ST  3pi 


+  _i  hV| £ 

[fi(L0)]2  3pi 
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Thus, 


3L° 

3L1 

3L° 

3L1 

-  (G_1) 

3F 

3H1 

3H1 

3L° 

1  h1^'  3L° 

Sq1 

3P,- 

3pi 

aq1 

n(L0) 

aq1 

api 

31 

3P, 

n  it 

!!(L6)  8»1. 

3F 

3H1 

3H1 

it!. 

1  h1^' 

aPi 

V 

31 

3q' 

n  n  u  t- 

n(L  )  3q  , 

(G-V  [  3F  3H1  _3F_  3H1  ' 

0(L°)  [W  3pi  ‘  3pi  Sq1  _  ’ 


But  (G*1)'  =  l/n(L°)  and  so  we  obtain  (3.1 ).■ 


The  cancellations  that  occur  to  yield  (3.1)  reflect  the  general  fact 
that  the  Poisson  brackets  before  and  after  reduction  correspond.  Similarly, 
if  K  is  a  function  of  (q,p)  we  obtain  the  formula 


(K,L'} 


5  tK-H'}’dr 


(K,F) 


(3.2) 


Thus,  if  K  is  a  first  integral  for  F,  then  (3.2)  becomes 


{K.L1}  =  -  1  {K.H1}  . 


(3.3) 


In  particular,  in  the  multidegree  of  freedom  case  in  which  all  but  the 
first  of  the  n  variables  (q,p)  are  in  action  angle  form  so  that 
(p^,  F}  =  0,  k  =  2,  ....  n  then  with  K  =  p^,  (3.2b)  becomes 

(Pk.  L1)--l(pk.H')=l^.  (3.4) 

3q 

Relations  of  this  type  are  useful  in  the  study  of  perturbations  of  inte- 
gs#ble' systems  when  n  >  2. 
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In  connection  with  the  identities  (3.1)-(3.4)  the  following  observa¬ 
tion  is  useful.  Along  an  orbit  for  the  unperturbed  system,  F  is  con¬ 
stant,  so  if  h  >  F,  will  not  vanish  and  — -U-  will  be  a  finite 

n(L°) 

constant.  Thus,  on  such  an  orbit,  (L  ,L^}  will  differ  only  by  a  multi¬ 
plicative  constant  from  {F,H^}. 

We  are  now  ready  to  state  our  main  result  for  Hamiltonian  perturba¬ 
tions  in  case  n  =  1. 


3 . 2  Theorem.  Consider  a  two  degree  of  freedom  Hamiltonian  system  of  the 
form 

He(q,p,-->,T)  «  F(q,p)  +  G(I)  +  eH1  (q.p,*,::)  ,  (3.5) 


and  assume  that  F  contains  a  homoclinic  orbit  (q(t  -  tQ),  p”(t  -  tQ)) 
connecting  a  hyperbolic  saddle  to  itself  (or  to  another  hyperbolic  saddle 
point).  Suppose  fl(I)  =  G'(I)  >  0  for  I  >  0.  Let  h-j  =  F(q',p')  be  the 
energy  of  the  homoclinic  orbit  and  let  h  >  h-j  and  =  G"\h  -  h^ )  be 
constants.  Let  (F,H^}(t  -  tg)  denote  the  Poisson  bracket  of  F(q,p)  and 

I  _  _ 

H  (q,p >.  evaluated  at  q(t  -  tQ)  and  p(t  -  tQ).  Define 


M(t0)  = 


x» 

(F,H1}(t  -  tQ)  dt ) 

’  -CO 


(3.6) 


and  assume  that  M(tQ)  has  simple  zeros  and  maxima  and  minima  of  0(1). 
Then  for  e  >  0  sufficiently  small  the  Hamiltonian  system  corresponding 
to  (3.5)  has  a  Smale  horseshoe  in  its  dynamics  on  the  energy  surface 
He  =  h. 

This  result  follows  from  our  previous  development  (the  reduction  and 
Proposition  3.1)  and  the  Melnikov  theory  given  in  the  references  at  the 


"S  B* 
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beginning  of  this  section.  Equation  (3.6)  can  also  be  obtained  from  the 
evolution  equation  for  F  along  the  unperturbed  orbit 

F  *  (F.H1 }  ,  (3.7) 

cf,  Arnold  [1964]. 

For  our  analysis  in  §5  we  shall  need  some  facts  about  the  construc- 
tructicn  of  the  horseshoe,  so  we  collect  them  here.  First  we  pick  an 

e  ®0  ®0 

energy  surface  H  =  h  >  h-j  and  consider  the  Poincard  map  P£  :  Z  Z 

(which  we  just  denote  P£  below)  associated  witn  the  periodically  per¬ 
turbed  system;  (2.9).  Here  I  =  { (q ,p , e )  1 0  =  9q  e  [0,2 tt]}  is  a  global 

cross  section  for  the  flow  of  (2.9).  By  hypothesis,  for  e  =  0  P£  has 

an  invariant  manifold  filled  with  a  continuous  family  of  (nontraverse) 
homoclinic  orbits.  If  M(tQ)  has  simple  zeros  then  this  manifold  breaks 
into  a  countable  set  of  homoclinic  orbits:  the  generic  case  found  in 
advanced  classical  mechanics  texts  (cf  Arnold  and  Avez  [1967],  Abraham 
and  Marsden  [1978]);  see  Figure  2a.  For  more  details  on  homoclinic  orbits 
of  maps  see  Moser  [1973]  or  Newhouse  [1980].  Here  we  merely  note  that 
the  Smale-Birkhoff  homoclinic  theorem  asserts  the  existence,  near  any 

transverse  homoclinic  point,  of  a  zero  difnensional  invariant  Cantor  set 

N 

A  or:  which  some  power  of  the  map,  P  ,  is  homeomorphic  tc  a  shift  on  two 

symbols.  Since  P^  |A  possesses  a  dense  orbit,  it  follows  (Moser  [1973]) 

that  (3.5)  possesses  no  analytic  second  integral. 

e* 

To  construct  the  horshoe  one  takes  a  small  rectangle,  R,  partially 
bounded  by  pieces  of  the  stable  and  unstable  manifolds  and  containing  a 
transverse  homoclinic  point,  integers^ ,  can  chosen  such  that  the 

forward  and  backward  images  P  '(R)  ,  P£  (R)  lie  in  a  neighborhood  U 


W 


of  the  saddle  point,  x.  The  linearized  map  DP£(x)  can  then  be  used  to 

approximate  the  motions  in  U  and  it  can  be  shown  that  there  are  an  integer 

N  <  ®  and  two  disjoint  'horizontal'  strips  C  P  (r)  =  b  whose  images 
N 

under  P£  are  disjoint  'vertical'  strips  V|  C  B  (Figure  2).  The  map 

PN:H.  -*■  V.  is  the  horseshoe, 
e  1  i 


(a)  The  perturbed  homoclinic 
orbit 


Wn(x) 


WS(x) 


(b)  Horizontal  and  vertical  strips 
Figure  2. 
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To  obtain  estimates  necessary  to  prove  hyperbol icity  of  A,  one  needs 
to  find  certain  sector  bundles  which  are  mapped  into  themselves  by  DP  .  In 
our  case  this  implies  that  the  choice  of  N  is  related  to  e  ,  the  pertur¬ 
bation  strength, si  nee  the  angle  between  the  tangent  vectors  of  the  manifold 
at  a  (transverse)  homoclinic  point  is  0(e)  (M(tg)  measures  the  0(e)  com¬ 
ponent  of  the  distance  between  the  perturbed  manifolds).  In  Appendix  B  we 
show  that  N  -  s,n(l/e).  Thus,  for  each  e  >  0  sufficiently  small  and  each 
h  >  h-j ,  there  is  an  invariant  set  near  every  transverse  homoclinic 
point  in  each  energy  surface  H£  =h  (cf.  Figure  4,  below).  However  as  e 
gets  smaller,  N  must  be  increased.  This  dependence  of  N  on  e  plays  an 
important  role  in  our  discussions  of  dissipative  perturbations  on  §5. 


4 .  Exampl e :  The  Coupled  Pendulum-Oscillator 


We  now  apply  Theorem  3.2  to  the  Hamiltonian  (1.4)  with  H1  given  by 
(1.10),  and  the  homoclinic  orbit  for  F  given  by  (1.7).  In  terms  of  the 
variables  (<J>,v,0,I),  we  have 


F(<M)  =  -  cos  4>  * 

and  H](<j),v,e,I)  =  ||  sin  0  -  ♦  . 

Thus 

{F.H1}  «  -v$  +  I-  v  sin  6  • 

(jl) 

The  energy  of  the  homoclinic  orbit  (1.7)  is  h^  *  1 ,  so  we  let  h  >  1 
let  =  ”  (h  -  !)•  Thus  (3.6)  gives 


M(tQ)  *  -  (4  sech(t  -  tg)arctan[sinh(t  -  tQ)] 


±  2l/2^"  —  sech(t  -  tQ)  sin(tui)}  dt  * 


The  first  term  is  odd  and  so  vanishes,  leaving 


co 

M(tQ)  =  ±  2>/2-^-  sech(t  -  tQ)  sin(toi)  dt  • 

•00 


This  is  evaluated  by  the  method  of  residues  as  in  Holmes  [1979],  yielding 


M(t0)  =  ±2 


Trail  . 

Tj  sin 


Since  M(tg)  has  simple  zeros  and  is  independent  of  c  we  conclude  that, 
for  e  >  0  sufficiently  small,  the  conditions  of  Theorem  3.2  are  satisfied 
and  we  have  horseshoes  in  the  Poincare  map  associated  with  the  pendulum- 
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oscillator  on  each  energy  surface  h  >  1,  where 

2 

He  =  tj-  -  cos  $  +  ojI  +  % 


sin  0  -  4> 


12 


=  h. 


Thus,  we  have  proved: 


4.1  Theorem.  The  Hamiltonian  system  with  energy  function 


2 

H(<j>,v,x,y)  =  -  cos  4>  +  j  (y2  +  w2x2)  +  |  (x  -  <(>)* 


has  horseshoes 
ciently  smal 1 , 


in  its  dynamics  on  each  energy  surface  H  >  1,  for 
and  hence  possesses  no  analytic  second  integral. 


suffi- 


J 
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5.  Non-Hamiltonian  Perturbations 


We  now  wish  to  consider  perturbations  under  which  the  total  energy  H 
is  not  conserved.  In  many  physical  problems  one  system  produces  energy  which 
is  subsequently  absorbed  by  a  second  system,  so  that  the  coupled  systems 
can  achieve  a  'dynamic  equilibrium'  in  which  (in  a  suitably  time-averaged 
sense)  energy  is  preserved.  This  often  manifests  itself  in  the  presence  of 
negative  damping  in  one  system  and  positive  damping  in  the  other.  We  will 
take  an  integrable  Hamiltonian  system  which  possesses  continuous  families  of 
non- transverse  homoclinic  orbits,  add  a  Hamiltonian  perturbation  He,  as 
before,  which  breaks  these  manifolds  to  give  transverse  homoclinic  orbits, 
and  then  add  dissipative  effects  which  cause  a  net  drift  in  the  energy  H 
for  perturbed  orbits  lying  near  the  homoclinic  manifold.  Under  suitable 
hypotheses,  such  dissipative  perturbations  can  leave  invariant  isolated  pieces 
of  the  continuous  family  Ah  of  horseshoes  discussed  in  Section  3.  For 
simplicity  we  shall  restrict  our  discussion  to  two  degree  of  freedom  systems 
(n  =  1). 

The  Hamiltonian  system  (1.9)  is  modified  to  include  dissipative  terms 
as  follows: 


•  9F  3H1  ,  f 

q  =  3^+  £lp"+  EVl 

3F  3H1  .  f 
P  =  '  aq  '  £  3—  +  £Y2f2 


6  =  J1(I)  +  e 


I  *  -e 


3H 

36 


+  e52g2 


(5.1) 


where  F,  Q  =  G'  and  H^  are  as  in  the  previous  sections  and  f^ ,  g^ 
are  functions  of  (q,p,I,0),  2ir  periodic  in  6.  Specific  hypotheses  on 
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f  ,  g.  and  on  the  dissipation  parameters  y^ ,  will  be  stated  subsequently. 

The  energy  function  H  =  H  -  F(q,p)  +  6(1)  +  eH  (q,p,I,9)  is  no  longer 
conserved,  and  our  earlier  equation  (2.1)  which  was  used  to  eliminate  I  is 
now  replaced  by 


H  =  e 


9He  .  .  3H£  .  ,  .  3He  „  .  .  3He  n  "1 

Y1  3q  f1  Y2  3p  f2  61  36  91  62  31  92J 


=  e 


3F 


3F 


Y1  3q  fl  +  y2  3p  f2  +  ^I)62g2j 


+  0(e‘) 


(5.2) 


=  eh(p,q,I)  +  0(e) 


Note  that  when  y..  *  5i  =0,  (5.2)  gives  H  =  0  and  (5.1)  becomes  (1.9)  with 
conservation  of  energy. 

The  five  equations  (5.1-2)  are  redundant  and  we  can  eliminate  the 
variable  I  by  regarding  I  as  a  function  of  p,  q,  0  and  H;  i.e.  by 
solving  H£  =  H  implicitly  for  I  but  remembering  that  H  is  a  variable 
with  its  own  evolution  equation  (5.2).  From  (2.2-5)  we  have 


I  =  L°(q,p,H)  +  eL1 (q,p,0,H)  +  0(e2)  (5.3) 

where 

L°  =  F_1 (H  -  F(q  ,p) ) ,  (5.4) 

and 

(q,P»8»l-  (q.P.H))  _  (5.5) 

fl(L°(q,p,H)) 


As  before,  we  have 

q'  =  q/9,  p'  =  p/0,  (5.6) 


and  from  the  implicit  equation  H£(q,p,0,L£)  =  H  we  obtain 


I 
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«fl 

v» 

*! 

i 

ii 


+  td  it!  =  o  etc 

ap  ai  3p  u*  etc>’ 


so  that 


3H”  .  - 

q  IF  "Vi 


becomes 


*  (if*  £Vi)/(lT  +  E5l9i) 

=  HIH  "(tH  vi])/(a(I)  4  e(It+  5i9i]) 


.  .  iki  .  Jitl .  Hi .  Mi)  +  ,,(  2, 

ap  (a  p  n  ap  a  ■'  K  }t 


A  similar  computation  for  p'  and  use  of  (5.2)  yields  the  three  dimensional 


system 


■  -  3L°  f Ikl.  Vi  3L0  6lgl]  O 

q  "  "  3p  1  3p  ■  fl  ‘  ap  a  J  +  0(e  )  , 


.  _  3L°  .  faL1  Y2f2  3L°  6lgl  n,  2, 
p  *  Iq"  +  e[1q  ‘  ~  aq~'lT  +  0(e  }  * 


(5.7a,b,c) 


01®  1  2 

H’  =efl[62g2  -  “Y^  -  %  Y2f2j  +  0(e  }* 


Equations  (5. 7a ,b,c)  constitute  the  system  we  now  study,  with  the  dependent 
variables  q,  p*  H  and  the  independent  time-like  variable  0.  For  y^  = 


6.  =  0  (5,7)  reduce  to  (2.9),  as  expected, 


To  deal  conveniently  with  the  slow  variable  H  compared  with  the 
fast  variables  (q,p)  we  use  a  slight  modification  of  the  usual  averaging 
theorem  in  which  the  0(6)  term  in  the  right  hand  side  of  (5.7c)  is  re¬ 


placed  by  its  8-average 


i . 


I 
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cfT  >  en(S2g2  -  ^Y,f,  -  %  Y2f2>: 


(5.8) 


1  2!tt 

where  an  overbar  denotes  the  time  average  ^  Iq  (  )  de.  The  modified 

averaging  procedure  is  discussed  in  Appendix  A.  Its  use  is  not  essential, 
but  it  makes  calculations  somewhat  easier.  Retaining  the  same  notation 


for  the  averaged  variables,  our  reduced  system  is  now  (5.7a,b)  plus 

H'  =  eh(q,p,H)  +  0(e3). 


(5.9) 


f-4 

A 


Let  6  =  denote  the  dissipation  coefficients  in  our 

®g 

system.  Also,  let  P  x  *  P_  .  denote  the  Poincar£  map  associated  With  the 

e,o  c»o  q 

system  (5.7a,b)  and  (5.9).  Thus,  P  maps  (an  open  subset  of)  IR3  to 

e  ,o 
3 

IR  and  is  given  by  advancing  the  independent  variable  9  by  2tt,  with 
starting  value  eQ. 

Let  us  assume  that  the  homoclinic  orbit  (q^e),  p\e))  for  F  has 

energy  =  h^  and  joins  a  hyperbolic  saddle  point  (q,  p")  to  itself.  (The 

case  of  a  heteroclinic  orbit  is  similar).  Thus,  for  each  value  of  H  > 

B0 

h1  ,  (q,p,H)  is  a  fixed  point  fcr  Pg  g*  Let  Q  denote  this  curve 

of  fixed  points.  Since  we  are  assuming  (q,p)  is  a  hyperbolic  fixed  point 

e0 

for  F,  Cq  q  is  a  hyperbolic  invariant  manifold  for  Pg  with  H 
restricted  to  an  interval,  say  h-|  <  Hg  <_  H  <_  .  Since  hyperbolic  manifolds 

are  preserved  under  perturbation,  we  have: 


5.1  Lemma.  For  <5  bounded  and  e  sufficiently  small,  there  are  invariant 


90  90 

curves  c  close  to  „  for  P  „.  Moreover,  the  stable  and  unstable 
-  e,S  -g -  0,0  - Q  e,6 - r- - 

manifolds  of  c  ® denoted  Ws(c  and  Wu(c  „)  are  "Cr  close"  to  those 

■  ”■  -  0  y  q  ~  £  f  0  _  £  J  O  — —  1  ” 


of  Cg  g  and  each  is  two  dimensional.  (See  Figure  3). 


Figure  3 


Notice  that  if  6=0  then  c  n  is  still  a  curve  of  fixed  points  by 

E  ,U 

conservation  of  energy  and  we  recover  the  situation  of  §3.  For  6/0, 

e0  a 

points  on  c  can  "drift"  under  iteration  of  P  %  since  energy  is  not 

e,6  q  e,6 

conserved.  However,  they  will  stay  on  c£  ^  (until  or  if  H  leaves  the  inter¬ 
val  [Hq.H,].) 

For  6=0,  suppose  Theorem  3.1  is  used  to  show  that  Ws(c  °n)  inter- 
sects  Wu(cr^n)  transversely.  This  persists  for  6  sufficiently  small,  by 
the  stability  of  transversal  intersections  under  perturbation.  Because  of 
potential  drift  in  the  H  variable,  this  alone  does  not  permit  us  to  con¬ 
clude  the  existence  of  horseshoes  for  6/0.  Rather,  we  must  control  H. 

The  crucial  hypothesis  that  enables  this  to  be  done  will  be  given  next. 


t'3i*  i 


( 
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Let  N  be  an  integer  fixed  so  that  (P£^q)M  has  a  horseshoe,  as 


0  \N 


described  following  Theorem  3.2.  Thus,  N  is  large  enough  so  that  (P  n) 
maps  2  horizontal  strips  in  B  back  around  to  two  vertical  strips  in  B 
again,  where  B  is  a  rectangle  lying  in  a  neighborhood  U  of  the  saddle- 
point  (Figure  2). 


I 


Let 


AH 


£ 


■ttN/S). 

h  dt  i 

-irM/fi 


(5.10) 


where  IT  is  given  by  (5.8)  and  h,  n  are  evaluated  on  the  homoclinic  orbit 
(q(9),  p(e))  at  an  energy  value  H.  (Recall  that  Q,  is  constant  on  this 
orbit).  From  (5.9)  we  see  that  AH  represents  the  approximate  change  in 
energy  in  following  a  point  starting  near  the  homoclinic  orbit  for  N  iter¬ 
ates;  i.e.  for  a  total  9-time  2ttN.  For  N  large  but  finite  P  .  maps 

&  9  0 

points  in  the  horizontal  strips  in  rectangle  B  in  Figure  2  back  to  the 
vertical  strips  in  rectangle  B  after  N  iterates.  Thus  AH  represents 
the  leading  term  in  the  energy  change  while  going  from  B  back  to  B.  Of 
course  AH  is  a  function  of  H  and  depends  on  e,  5  and  N.  Strictly 
speaking,  (5.10)  should  be  evaluated  on  trajectories  just  inside  the  homo¬ 
clinic  orbit,  but  as  we  show  in  Appendix  B,  as  e  0,  N  -+•  ®  and  the  hori¬ 
zontal  and  vertical  strips  H.,  V.  c  b  must  be  taken  closer  and  closer  to  the 
homoclinic  orbit.  Thus,  since  we  only  need  AH  to  leading  order  in  what 
follows,  evaluation  on  the  homoclinic  orbit  is  sufficient. 

Now  we  state  our  basic  energy-transfer  condition: 

Condition  (H).  Assume  there  is  a  value  Hc  >  1  of  H  at  which  AH  given 
by  (5.10)  changes  sign  transversely,  i.e. 
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and 


AH  <  0 

11 

H  <  H 

c 

AH  >  0 

if 

H  >  H  , 

c 

Ih. 

/  0. 

Under  this  condition  and  e  sufficiently  small,  we  have 


5.2  Lemma.  There  is  a  smooth  function  H(q,p)  defined  on  the  rectangle  B 

0 

such  that  if  (q , p )  €  B  then  (P  ® JN(q,p,H(q,p) )  has  the  form  (q,p, 

c  )  O 

H(q,p)). 


O 

Proof.  Since  the  exact  energy  change  differs  from  (5.11)  by  0(e),  persis¬ 
tence  of  transversal ity  guarantees  that  condition  (H)  is  also  true  for 
the  exact  energy  change  for  e  sufficiently  small.  For  6=0,  the  surface 
H  =  constant  =  H  is  preserved  by  (P  un)  .  For  6^0  the  surface  W 

C  O  y  U  V 

is  preserved  and  contracting  to  first  order  in  e.  By  persistence  of  hyper¬ 
bolic  invariant  manifolds,  there  is  a  nearby  surface  exactly  invariant  under 
(P^0^)'  ;  this  surface  is  the  graph  of  H.* 

Thus,  we  have  identified  a  surface,  say  I  near  H  =  H  such  that 
0  c  c 

(Pe“/  maps  Ec  to  £c. 

Now  we  wish  to  show  that  there  is  a  horseshoe  in  this  surface  zc.  By 
the  arguments  in  Holmes  and  Marsden  [1981,  Appendix  A],  we  must  check  that 

Ws(c^0  )  and  Wu(cGe>  )  continue  to  intersect  transversally,  for  6^0. 

£  j  (5  £  >6 

To  do  this,  we  form  the  Melnikov  function  at  energy  value  Hc  for  the 
system  (5.7)  and  (5.8). 

If  M(tQ)  is  given  by  (3.6),  then,  using  Proposition  3.1,  the  Melnikov 
function  for  (5.7a)  and  (5.7b)  is  given  by 
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1 

f  3L°  f 

SL°  f  } 

M  a 

'  —OO 

l  3  p  Y2f2 

(Note  that  the  6^2  term  cance1s  out).  Thus,  using  (5.4),  we  get 


5.3  Lemma: 


W  =  -2  M^o)  + 


3F  ,  3F  ,  ' 
3p  Ylf1  '  3q  Y2f2 


•] 


(t  -  tn)  dt  (5.11) 


We  will  assume  that  6  is  chosen  such  that  M  (tg)  continues  to  have  simple 
zeros. 

In  Appendix  B  we  discuss  the  relationship  between  N  and  e.  As  e 
gets  smaller,  the  number  of  iterates  required  of  the  Poincare  map  to  guarantee 
a  horseshoe  gets  larger.  It  is  shown  that: 


5.4  Lemma:  We  have 


N  =  L 


+ 


2a  In 


B 


7 


(5.12) 


where  a  and  B  are  constants,  L  is  a  fixed  integer  and  M(Hc)  is  the 
supremum  of  M(tg)  over  tg;  (the  dependence  of  M  on  Hc  is  made  explicit 
in  (5.12)). 

This  result  applies  to  the  case  of  Hamiltonian  perturbations  (5  =  0).  When 
5^0  there  isan  analogous  result  N  =  N(e,6)  in  which  a  =  a(6)  and  M(HC) 
is  replaced  by  Mg(Hc).  However,  in  our  application  we  set  6  =  to  be 
0(eM),  so  that  e6  =  <<  e  for  e  <<  1  and  the  effects  of  6  in  (5.12) 

can  be  ignored,  cf.  Appendix  B. 

From  (5.11  )we  obtain  a  condition  on  the  size  of  y^  and  Y£i  if  M(tQ) 
has  simple  zeros  and  oscillates  with  an  amplitude  M(HC),  then  (5.11)  gives 
us  conditions  of  the  form 
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c]Yl  <  M(Hc) 
c2Y2  <  M(HC) 

which  must  be  satisfied  for  M5  to  still  have  simple  zeros. 
We  summarize  our  findings  as  follows: 


(5.13) 


5.5  Theorem.  Suppose  that  N,  H„  and  5  can  be  chosen  so  that  conditions 

^  0 

(H),  (5.12)  and  (5.13)  all  hold.  Then  some  iterate  (P  of  the  Poincarg 

map  of  the  reduced  system  (5.7)  has,  for  e  sufficiently  small,  a  horseshoe 
in  its  dynamics;  the  horseshoe  lies  near  the  homoclinic  orbit  in  the  (q,p) 
variables  and  near  the  (non- invariant)  energy  surface  H  =  Hc. 
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6.  Example:  The  Oscillator-Pendulum  with  Positive  and  Negative  Damping 


We  now  wish  to  show  that  for  e  sufficiently  small  the  pendulum-oscilla¬ 
tor  system  considered  in  §4, 


4>  =  v,  v  =  -sin  +  e(x  -  4>)  * 
v  =  y,  y  =  -O)  x  +  e(4>  -  x), 

continues  to  have  a  horseshoe  when  dissipation  is  included.  Specifically, 
we  add  negative  damping  (-6)  to  the  oscillator  so  that  it  drives  the  pendulum, 
which  now  has  positive  damping  (y): 


4>  =  v  , 

v  =  -sin  4>  +  e(x  -  4>)  -  cyv  , 
x  =  y  , 

y  =  -uj2x  +  e(4>  -  x)  +  e6y  , 


In  action-angle  variables,  (6.1)  becomes 


4>  =  v  . 

v  =  -sin  <p  +  e 


/ 21  . 

/—  sin  0  -  <p 

'  U) 


-  cyv  , 


9  =  w  +  e 


f  /2f  . 

1  /—  sin  0  -  <j> 
/  uj 


]  sin  9 


eg  sin  9  cos  0  . 


I  =  -e 


/Ssin 

V  OJ 


e  -  4) 


JITu 

/2Ito  cos  0  +  &S2I  cos20  .  ^ 


(6.1) 


(6.2) 


Note  that  6  >  0  represents  damping  while  6  >  0  represents  negative  damping 
(energy  production).  The  energy  evolution  equation  (5.2)  is 


^  =  efcy2  -  yv2)  , 

H'  =  (26u)I  cos20  -  yv2)  +  0(e2).  (6.3) 


28 


where  the  dependent  variable  is  now  6.  Using 

I  =  L°(<M,H)  =  l  (H  -  (v2/2  -  cos  4>)),  (6.4) 

from  (2.4),  (6.3)  becomes 

H'  =  |  [6(H  +  cos  ct>  -  v2/2) (1  +  cos  20)  -  yv2]  +  0(e2).  (6.5) 

Although  we  do  not  need  them  explicitly  in  the  calculations  to  follow,  we  also 
give  the  reduced  evolution  equations  for  <j> ,  v: 


1  L  +  sfy [<Lsin_e  _  sin2 9 

w  L  “[  K2A  u 

1  F-  sin  *  +  £  (sin  *  fsla 


-  -  -1  ■— 9  +  6  sin  6  cos  e) )1  +  0(e2)  7 

O)  \  '  1 


■0  4>  sin  0 


y  sin  0  cos  0  + 


+  (/2 sin  e  -  o>4>  -  yv)j  +  0(e2),  (6.6a,b) 


where  A  =  H  +  cos  4>  -  v2/2.  Equations  (6.5a,b)  and  (6.5)  correspond  to 
(5.7a ,b,c) . 

We  now  average  (6.5).  The  transformation  (4>,v,H)  ■>  (4>,v,F)  is  given 
by  (4>,v,H)  =  (4>,v,H  +  cu(4>,v,F,0)),  where 


(6.9) 
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We  now  check  condition  (H).  Inserting  the  unperturbed  solution  (1.7)  we 

2, 


have  H  +  cos  <J>(t)  -  v(t)4/2  =  H  -  1  and  (6.9)  becomes 

H'  =  J  r«(H  -  1)  -  4y  sech2|]  +  0(e2). 

H  -  e[6(H  -  1)  -  4y  sech2t]  +  0(e2). 


or 


(6.10) 


Hence 


7tN/u) 


AH  =  el  [<S(H  -  1)  -  4y  secfrt]  dt  . 
1  -irN/w  7 


.  E[m  (H  .  D  .  8V  ta„h(£)]  , 


(6.11) 


and  so  condition  (H)  is  satisfied  if  we  have 


H  _  i  +  8<jj  tanh (ttN/oj) 
'c  2ttN 


X 
( 6) 


Sk  6 (- or  N  large) 


(6.12) 


We  next  compute  the  Melnikov  function  ^(tQ)  from  (5.11). 

Section  4  and  (3.1)  we  have 


From 


M(tn)  =t\  2tt/2(H  -  1)  sech 
U  <0^ 


r 

TT(0 
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-oH  sin  to  t„  . 


(6.13) 


O 

Using  F  *  v  /2  -  cos  <p »  Y i  f i  =  0»  and  ^^2  =  10  (5.11),  the  second 

term  of  M^jY)(t0)  is 


.00 


_oo  U)“ 


-r  v(-*v)  dt  =  -J  (2  sech(t  *tn))2  dt  =  ^  . 

J  aa  U  .  .  C 


U) 


Thus  we  obtain 


(S,Y)(tO>  -  1 )  seoh(^.j  sin  to  tg  +  8y] 


M 


(0 


def  M(H)  sin  w  tg  +  8y/u> . 


(6.14) 


To  compute  our  verification  of  the  hypotheses  of  Theorem  5.5,  we  note 
that,  for  M^tg)  to  have  simple  zeros,  we  require 


ll 


I 
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M(Hc)>8y/w% 

- —  V2: 

tt  sec!r(W2)  ' 


while  for  satisfaction  of  condition  (H)  we  have 


u  -  i  =  Jux.  def  c  JL.  . 
nc  1  2ttN<5  -  u2  N6 


(6.15) 


(6.16) 


We  also  have  the  relationship 


i 


N  =  N(e)  =  L  +2aAn(S/e/Hc  -  1) 


(6.17) 


from  (5.12),  where  the  remaining  constants  in  ^(tQ)  are  accumulate<l  1nto 
0.  From  (6.16)  and  (6.17)  we  have 


(Hc  -  1)[  L  +  2a  Jln(B/e/Hc  -  1)]  *  0^/6 


(6.18) 


To  satisfy  (6.15)  and  (6.18)  simultaneously  we  pick  y  and  6  to  be  of 
the  same  order  in  e,  say  )=  {  5  cp{;  y  >  0.  Then  (6.18)  may  be 


rewritten  as 


(Hc  -  1)[1  +  C3  iln(6/c/Hc  -  1)]  =  C4 


(6.19) 


A  simple  exercise  in  calculus  shows  that, for  e  small  (6.19)  has  a  unique 


solution  Hc  near  Hc  =  1  and  for  any  fixed  a,  0  <  a  <  2, 


H  -  1  >  ea, 


for  small  e;  "how  small"  depends  on  how  close  °  is  to  0. 


(Proof.  Rewrite  (6.19)  as 


(6.20) 


»  ••  < 
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e^f - 


f 


where  x  =  H  -  1.  By  considering  the  graphs  of  y  =  ev/x  and  y  =  <p (x ) 
one  sees  that  for  small  e,  (6.19)  has  a  unique  small  solution  x(e)  0 
as  e  ->  0.  Moreover,  as  $  vanishes  to  all  orders  at  x  =  0,  <j>(x)  <_  x15 
for  p  >_  1/2  and  x  small.  It  follows  that  x(e)  is  larger  than  the 
solution  of  e/x  =  x*3  for  small  £‘,  i.e.  x(e)  >_  £2/(2p-l)  for  sman 
Let  a  =  2/(2p-l ).  ■) 

Picking  c  <  2y,  (6.15)  is  now  easily  satisfied,  since 

H  -  1  >  c  e°  >  Cie^y^  (6.21) 

cal 


for  e  sufficiently  small.  For  example,  we  can  take  y  =  1/2,  since  then 
the  damping  perturbations  e(ey6),  e(eVy)  appear  at  0(e3//2)  and  the  O(e^) 
terms  ignored  in  our  computations  do  not  affect  the  results.  Thus,  we  have 
proved: 


6.1  Theorem.  The  system  (6.1)  has  a  horseshoe  in  its  dynamics  provided 
3/2—  3/2— 

we  choose  ey  =  e  y  and  e5  =  e '  6  and  e  sufficiently  small. 

Of  course  it  is  possible  to  vary  the  orders  of  6  and  y  with  some 
latitude  and  still  maintain  the  hypotheses.  Specifically,  6.1  remains  valid 
if  we  choose  y  =  evy  and  6  =  ,0<_y£v<  1,  v  j*  0. 
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Conclusions 

In  this  paper  we  have  developed  applicable  techniques  for  establishing 
the  existence  of  chaotic  dynamics  in  the  sense  of  the  presence  of  a  horse¬ 
shoe  for  both  Hamiltonian  and  non-Hamiltonian  perturbations  of  systems  with 
two  degrees  of  freedom  containing  homoclinic  orbits  and  periodic  orbits. 

While  horseshoes  are  not  strange  attractors,  they  are  often  visible  and 
behave  like  them  in  numerical  experiments  (perhaps  due  to  small  background 
noise);  cf.  Franks  [1981].  Arnold  diffusion  is  a  higher  dimensional  manifes¬ 
tation  of  the  same  phenomenon  and  is  certainly  seen  in  many  examples  (see, 
for  instance  Lieberman  [1980]). 

For  conservative  perturbations  the  method  is  a  straightforward  combina¬ 
tion  of  a  classical  reduction  scheme  with  a  method  of  Melnikov.  For  non¬ 
conservative  perturbations  a  delicate  energy  balance  argument  is  needed  to 
ensure  that  at  least  one  horseshoe  survives  near  the  energy  balance  point. 

Near  other  points  there  is  a  "ghost  horseshoe"  which  decays  because  of  energy 
drift.  If  the  dissipation  terms  all  contribute  to  energy  loss  then,  while 
no  invariant  set  remains  near  the  homoclinic  orbit  (since  H  decreases  on 
all  orbits),  the  manifolds  Ws(c^c  J  wu(c®rt)  continue  to  intersect  and 

£  »  0  £  9  0 

the  resulting  ghost  horseshoes  would  give  rise  to  complicated  dynamics  on 
finite  time  intervals,  as  orbits  move  through  the  energy  band. 

The  results  are  shown  to  apply  to  typical  perturbations  of  the  pendulum- 
oscillator  system,  thereby  showing  that  this  classical  example  has  complex 
dynamics  and,  in  particular,  is  non-integrable. 
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Appendix  A.  A  Modified  Averaging  Theorem 

For  the  basic  averaging  theorem  see  Hale  [1969].  Here  we  consider  a 
system  of  the  form 


x  =  f(x)  +  eg(x,y,t), 
y  =  eh(x,y,t),  0  <  e  «  1 , 


(A.  1 ) 


where  x  =  x(t)  and  y  =  y(t)  are  the  fast  and  slow  variables,  the  functions 
f,  g,  h  are  sufficiently  smooth  and  the  latter  two  are  T-periodic  in  t. 


A1 .  Proposition.  There  exists  a  near  identity  time  dependent  change  of  co¬ 
ordinates  (x,y)  (x,z)  under  which  (A.l)  becomes 


x  =  f(x)  +  eg(x,z,t)  +  0(e2) 
z  «  eF(x,z)  +  0(e2)  , 


where  h 


h(x,y,t)  dt 


is  the  t-average  of  h. 


(A. 2) 


Proof.  As  in  the  usual  averaging  theorem,  we  set 

y  =  z  +  eu(x,z,t)  .  (A. 3) 

Differentiating  (A. 3)  with  respect  to  time,  we  obtain 


y  =  z  +  eu  +  eD  ux  +  eD  uz  ,  where  ()  =  3/3t.  (A. 4) 

X  i 

Using  (A.l),  (A. 2)  and  (A. 4): 


(Id  -  eD2u)z  =  y  -  eu  -  eDxtix  =  e[h(x,z  +  eu,t) 
-  eu(x,z  +  eu,t)  -  Dxuf(x)] 


(A. 5) 


3H 


We  can  write 

h(x,z,t)  =  h(x,z)  +  ii(x,z,t)  ,  (A. 6) 

where  h  is  T-periodic  in  t  and  for  zero  mean.  From  (A. 5-6)  we  have 

(Id  -  eOzu)z  =  e[h(x,z)  +  h(x,z,t)  -  u(x,z,t)  -  Dxu(x,z,t)f (x)] 


+  0(e2)  • 

(A. 7) 

Thus,  if  we  set 

§£  +  Dxuf(x)  =  h(x,z,t)  , 

(A. 8) 

we  have,  from  (A. 7) 

z  =  eh(x,z)  +  0(e2)  ? 

(A. 9) 

and,  using  (A.2)  in 

(A.l) 

x  =  f (x)  +  cg(x,z,t)  +  0(e2) 

(A. 10) 

It  remains  only  to  check  that  the  linear  partial  differential  equation 
(A. 8)  admits  a  solution  u  =  u(x,z,t).  However,  (A. 8)  has  the  solution 

(t  , 

u(x(t),z,t)  =  u(x,z,0)  +  h(x(s) ,z,s)  ds  (A. 11) 

Jo 

where  x(t)  satisfies  x(0)  =  x  and  x  =  f(x).  ■ 

Most  of  the  usual  averaging  results  go  through;  in  particular,  solu¬ 
tions  (x(t),  z ( t ) )  of  (A.2)  remain  v/ithin  0(£2)  of  those  of  (A.l)  for  times 
of  0(l/e).  Since  we  wish  only  to  integrate  for  times  of  0 ( N )  =  0(£n(l/e)) 
(eqs.  (5.10)(5.12))  the  averaged  equation  may  be  used  in  computations.  Mote 
that  the  transformation  u(x,z,t)  is  not  in  general  T-periodic,  since  solu¬ 
tions  of  (A. 8)  depend  upon  the  (nonperiodic)  flow  x(t)  of  x  =  f(x).  One 
must  therefore  be  careful  in  inferring  the  existence  of  T-periodic  solutions 
of  (A.l)  corresponding  to  fixed  points  of  (A.2),  as  in  the  usual  averaging 
theorem. 
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Appendix  B.  The  Iteration  Number:  Proof  of  Lemma  5.4 

In  this  appendix  we  derive  a  relationship  between  N,  the  number  of 

e  o  0  n  m 

iterates  of  the  Poincare  map  P  n  necessary  to  guarantee  that  F  =  (P  u  )r 

has  a  horseshoe,  and  s,  the  perturbation  parameter.  Henceforth  we  drop 

®0 

the  sub-  and  superscripts  on  P  n. 

Let  x  be  the  (perturbed)  saddle  point  of  P  and  y  a  transverse 
homoclinic  point  lying  outside  a  ball  B^(x)  of  radius  p  about  x.  The 
Melnikov  theory  tells  us  that  the  maximum  distance  betv/een  the  manifolds  near 
y  is 

dmax  =  eK1H<Hc>  +  I6'1) 

where  M(H  )  =  sup  M(tn)  and  K-,  is  a  constant. 

C  tQS[0,T)  0 

We  next  need  a  basic  result  from  dynamical  systems  theory,  the  "lambda- 
lemma1'  (Palis  [1969],  Newhouse  [1980]),  which  enables  us  to  make  our  choices 
of  horizontal  and  vertical  strips  in  the  horseshoe  map  more  precise: 


B.l  Lemma.  Let  x  be  a  hyperbolic  saddle  point  of  a  Cr 
P  and  D  c  Wu(x)  an  open  disc  in  its  unstable  manifold, 
of  dim(Wu(x))  meeting  Ws(x)  transversely  at  a  point  y. 
contains  discs  arbitrarily  Cr  close  to  Dy. 


diffeomorphism 

Let  A  be  a  disc 

Then  U  Pn(A) 
n  >  0 


This  result  implies  that,  if  y  e  Wu(x)  Ws(x)  is  a  transverse  homo¬ 
clinic  point,  then  Wu(x)  and  Ws(x)  accumulate  on  themselves,  giving  us 
the  structure  of  Figure  B.l.  We  assume  that  the  map  is  orientation  preser¬ 
ving,  as  are  the  Poincar4  maps  occurring  in  the  application  of  this  paper. 
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a 


?* 


We  pick  a  rectangle  R  bounded  by  pieces  of  Ws(x),  Wu(x)  as  shown. 

Since  y  is^fixed  independent  of  e,  there  are  fixed  integers  L] ,  L 2 
such  that  P  ^(R),  P  ^(R)  C  B^(x)  and  the  'height'  and  'width'  of 
P  (R)  and  P  ^(R)  are  eK»M(H  ),  eK4M(H  )  respectively.  Once  in  B  the 
dynamics  is  dominated  by  the  linearized  map,  which,  working  in  suitable  co¬ 
ordinates,  we  can  take  to  be 

np  =  °)  u,  i  ^  i  ^  t,i 


DP  -  A  J  .  M  <  1  <  |A 

0  Y 


-(Vo+N«)  L.+N. 

To  obtain  the  horseshoe  structure  P  c  c  (R)  n  p  1  '(R)  as  shown,  we 

require  further  iteration  numbers  M.  =  N.(e)  such  that 


"l 

JJ  eK3M(Hc>  ■  ^  > 


£K4M(hc)  =  Kgy 


(e)  =  a«,n(gi/eM(Hc))  i-1,2, 
a  =  l/£n(y/X  ) . 


Thus  the  total  number  of  iterates  is  N  =  L-j  +  Lg  +  Ni  +  or 


N  =  N(e)  =  L  +  2a£n(g/EM(Hc)). 


where  L  is  a  fixed  integer  and  a,  B  are  constants. 

When  5^0  M(HC)  should  be  replaced  by  M5(HC)  and  X,  y  by 
^  +  Ky6 1  y  +  Kg5,  leading  to 

^-(M)  =  a(6)£n(Bi/cM6(Hc))  , 


(B.2) 


(B.3) 


(B.4) 


where 


and  Kg  is  a  positive  constant.  Thus  N(e,6)  >  N(e)  in  general.  However, 

i  to 

in  our  application  we  take  6  of  order  cy,  y  >  0  (for  example  e '  )  and 

thus  e6  =  e1+u6,  say,  and  the  dependence  of  N  on  6  is  weaker  than  its 
dependence  on  e,  and  hence  can  effectively  be  ignored  in  thelimit  e  -*•  0. 
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